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ABSTRACT 

We use lubrication theory and matched asymptotic expansions to model the quasi-steady 
propagation of a liquid bridge through an elastic tube. In the limit of small capillary number, asymptotic 
expressions are found for the pressure drop across the bridge and the thickness of the liquid film left 
behind, as functions of the capillary number, the thickness of the liquid lining ahead of the bridge and the 
elastic characteristics of the tube wall. For a given precursor thickness, we find a critical propagation 
speed, and hence a critical imposed pressure drop, above which the bridge will eventually burst, and 
hence the tube will reopen. 


INTRODUCTION 

An airway in a lung may crudely be described as an elastic tube with a viscous liquid lining. It is 
well-known that this lining is susceptible to an instability driven by capillarity at its free surface. If the 
initial lining thickness is sufficiently large, the instability culminates in the formation of a liquid bridge 
(known in the medical context as a bolus) which blocks the tube. This phenomenon, known as airway 
closure, can lead to respiratory difficulty. It is known to be exacerbated by increased liquid content in the 
lung and by increased airway wall flexibility (ref. 1 ). In normal gravity, airway closure is generally 
confined to the lower regions of the lung. However, in microgravity airway closure still occurs and 
appears to be more homogeneously distributed (ref. 2). 

Once such a liquid bolus has been formed, it propagates along the airway under the pressure 
drops imposed across it during respiration. As it does so, it leaves behind a thin liquid film. If the 
thickness of the trailing film exceeds that of the liquid lining ahead of the bolus, then the bolus must 
decrease in volume as it propagates, leading ultimately to reopening of the airway. In this paper we model 
the propagation of a liquid bolus along an elastic tube and hence deduce minimum criteria for airway 
reopening. In particular, we find that a more compliant airway is easier to reopen. 

SYMBOLS 


See figure 1 for a definition sketch. 
a tube radius 

U bolus propagation speed 

o liquid surface tension 

p liquid viscosity 

h h 2 trailing film thickness, precursor film thickness 

P h P 2 , P h pressure behind bolus, ahead of bolus and in bolus respectively. 

w tube wall displacement 

Ed tube wall stiffness (Youngs modulus times thickness) 
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Dimensionless parameters 


Capillary number 
Dimensionless pressure drop 

Dimensionless wall flexibility 


Ca = \x Ulo 


A P 
G 


a (P r P 2 ) / 2 o 


o 

Yd 



(3 Caf 3 


ASYMPTOTIC ANALYSIS 


Throughout this paper we assume that the flow is axisymmetric and that inertia may be neglected 
(i.e. that the Reynolds number is small); our governing equations are therefore the axisymmetric Stokes 
equations. At the free surface of the liquid, we have the classical balance between liquid stress and inter- 
facial tension. We perform an asymptotic analysis of the problem, assuming that the capillary number is a 
small parameter. This approach has been adopted by many previous authors, and applied to such 
problems as the flow of bubbles in tubes (ref. 3, ref. 4), the dynamics of free surfaces in Hele-Shaw cells 
(ref. 5) and plate withdrawal (ref. 6). Our analysis represents a generalisation of these papers to allow for 
flexibility of the tube wall. In this paper we employ a very simple (essentially linear) relation between the 
inwards wall displacement w and the pressure drop across it, namely 


Py 



Ed w 
(a-w) a 


(i) 


The key observation is that at small capillary number, in most of the flow domain viscous 
dissipation is dominated by surface tension. Indeed in the leading-order problem, obtained by setting the 
capillary number to zero, there is no flow and the problem is simply one of capillary statics. Formally, 
nondimensionalising the Stokes equations (using (7, a and 2o/a as scales for velocity, length and pressure 
respectively) reveals that 


P b ~ const. + Q(Ca\ 


( 2 ) 


and so the two free surfaces are, to leading order, surfaces of constant mean curvature, that curvature 
being the ratio between the constant pressure difference across each surface and the interfacial tension o. 
Also imposing that the surfaces be axisymmetric and analytic, the only admissible constant-mean- 
curvature surfaces are hemispheres, with leading-order mean curvature 2/a. 

We can anticipate that the scalings employed above will break down where the liquid free 
surfaces become close to the tube wall, since the hemispheres described above cannot be joined 
analytically to the uniform films ahead of and behind the bolus. Instead, there are transition regions (see 
figure 1 ) between the uniform films and the hemispheres, in which the pressure ceases to be constant as 
viscous dissipation becomes important. However, asymptotic simplification can still be achieved since in 
the transition regions the liquid layer is thin, with a slowly-varying free surface, so that we can employ 
lubrication theory. 

Consider first the rear meniscus. If the film thickness h and wall displacement w are 
nondimensionalised with h, and axial distance z with /?, (3Ca) J \ then in a frame moving with the bolus h 
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( 3 ) 


is found at leading order to satisfy a variant of the Landau-Levich equation: 

w"'(z) * h'"(z) = — . 

h 3 


Coupled to this is the leading-order dimensionless wall law from (1), namely 

w = G(w " + h"). 


( 4 ) 


where G is the dimensionless flexibility: as G - 0, the tube becomes rigid and the classical Landau- 
Levich equation is recovered. 

Equations (3) and (4) combine to a single ordinary differential equation for q := h + w (i.e. the 
free-surface profile): 


»1 


m 


q - Gq" - 1 
(q - Gq") 3 ' 


( 5 ) 


An initial condition for (5) is that the film become uniform far away from the bolus, and by linearising 
about q = 1 it is apparent that this one initial condition specifies the problem completely. Moreover the 
unique solution satisfying this condition has the asymptotic behaviour 

q - 1 + e lz as z - - <*>, (6) 

where / is the unique real, positive root of the cubic 

X 3 + Gk 2 - 1 = 0; (7) 


equation (6) is used as an initial condition in integrating (5) numerically. 

The solution of (5) with initial condition (6) is found to behave quadratically for large, positive z: 


n 


1 7 

— Az + Bz + C as z 

2 


( 8 ) 


where A , B and C are numerically-determined constants. Notice that these constants are not all uniquely 
determined because the origin forz may be chosen arbitrarily; however /I and the combination 
F = AC-B 2 / 2 are independent of the origin chosen for z and so are determined uniquely for any fixed G. 

The final step is to match the transition region with the outer, bolus solution. Formally, we apply 
Van Dyke’s matching rule (ref. 7) to the three-term outer solution (in which the free surface is 
hemispherical up to O(Ca) ) and the one-term inner solution found above. Matching the curvature gives 
the well-known asymptotic relation for the trailing film thickness: 

— - A(G)OCaf 3 . (9) 

a 
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Moreover, the higher-order matching gives asymptotic expressions for the apparent contact radius and 
contact angle of the meniscus. From these, the first perturbation to the meniscus curvature can be found, 
and hence the pressure difference across the rear meniscus is given asymptotically by 

~ p „) ~ 1 + F(G)QCa)™. (10) 


Now we apply the same asymptotic arguments to the front meniscus. The same differential 
equation (5) is found for (here nondimensionalised with the precursor film thickness h 2 ), with G 
replaced by G := G h x /h 2 and the initial condition of uniform r| now imposed for large positive z. In 
contrast to the rear meniscus, this problem admits a one-parameter family of solutions, with asymptotic 
behaviour 


T) ~ 1 + a e ' zj 12,-1 cos (\li4l 2 - 1 )z/2l' 2 ) as z - *>, 


( 11 ) 


(/is the real positive solution of (7) with G replaced by G). Once again the transition film behaves 
quadratically as it approaches the meniscus: 

1 ~ ~ 

n - — Az +Bz + C as z 
2 


( 12 ) 


Now, however, the translation-invariant constants A and F = AC -B 2 tlare both functions of the initial 
condition a. Matching the parabola (12) with the front meniscus gives 


A 


a(3Ca) 


2/3 


( 13 ) 


the shooting parameter a must be adjusted until A satisfies this (recall that h 2 is a physically specified 
quantity, while A, is a quantity to be found). The pressure difference across the front meniscus is given 
by 


T-(^2 - P o) ~ 1 + F (3 Ca) m . 


( 14 ) 


Here, F is in effect a function of A : for each fixed G they may be plotted against one another using a as 
a parameter. Hence the total pressure drop driving the bolus is related to its propagation speed by 


a 

2a 


(*\ 


" *a> 



- h 2 )} 

F(G) - F 

G ; 2 


a{3Cd) m ) 


(3 Ca) 


in 


( 15 ) 


For a rigid tube, when G = G = 0, this means that, given the precursor layer thickness h 2 and 
the imposed pressure drop P } -P 2 , the propagation speed of the bolus is given by (15); then the trailing 
film thickness h ] can be found from (9). In general the situation is more complicated; (15) and (9) are 
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coupled since /i, appears in G. 


REOPENING CRITERION 


The condition for the bolus eventually to rupture is that h { be greater than h 2 . From (9) we can 
deduce a critical capillary number above which this condition is satisfied: 


reopening 


Ca > Ca c 


h ) 

3/2 

{ 

\ 

"2 

/ 

a 

o 



, a , 

[h,Edj 


( 16 ) 


where/is a numerically-determined function. Hence from (15) the minimum pressure difference A P c 
required to reopen the tube takes the form 


A P 


C 


( a o 

~^ g \h 2 Jdy 


( 17 ) 


where g is a second numerically-determined function. 

Typical critical pressure drops are plotted against the flexibility parameter a/Ed in figure 2. It 
can be seen that increased wall flexibility decreases the critical pressure drop, i.e. makes the tube easier to 
reopen. The effect appears to be very slight compared to the increase in critical pressure drop caused by 
increasing precursor layer thickness; this is not surprising given that the present analysis is only valid for 
relatively stiff walls. 
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Figure 1 : Definition sketch of a bolus in a flexible tube. 



Figure 2: Critical dimensionless pressure drop (A P) c = a(P l -P 2 )/lo versus flexibility 
parameter o/Ed for various values of the dimensionless precursor thickness. 
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Convective Flows 




